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Abstract. On January 19, 1993, a very bright bolide 
(peak magnitude —23) crossed the sky of Northern Italy, 
ending with an explosion approximately over the town of 
Lugo (Emilia Romagna, Italy). The explosion (14 kton 
of energy) generated shock waves which were recorded by 
six local seismic stations. A reanalysis of the available data 
leads us to the hypothesis that the meteoroid was a porous 
carbonaceous chondrite, somehow similar in constitution 
to the asteroid 253 Mathilde. 

Key words: meteors, meteoroids - minor planets 



1. Introduction 

The atmospheric interaction of large meteoroids provides 
our primary tool to characterize their population, physi- 
cal and chemical properties, and dynamical evolution. In 
turn, this can lead to a better understanding of the diverse 
populations of small bodies of the Solar System. Currently, 
our knowledge is still quite limited, although, especially af- 
ter the impact of comet D/ Shoemaker-Levy 9 on Jupiter, 
the research efforts in this field have been intensified. In 
particular, in 1994 the US Department of Defense made 
of public domain its records on energetic bolides over a 



time span of about twenty years (Tagliaferri et al. 1994). 
These data indicate that, from 1975 to 1992, there were 
136 airbursts of energy greater than 1 kton, but the real 
number was probably at least 10 times higher, because the 
satellite system does not cover the entire Earth surface. 

Both data and theories are required to assess the im- 
pact hazard and to understand the very bright bolides. 
From this point of view, the Lugo bolide is a very inter- 
esting event, because the airburst was detected by several 
seismic stations. The corresponding data allow us to char- 
acterize the meteoroid and to draw some tentative infer- 
ences about its nature and origin. We have carried out a 
reanalysis of this event and we found that the data are 
most consistent with the hypothesis that the involved me- 



teoroid was a porous carbonaceous chondrite, somehow 
similar to the asteroid 253 Mathilde. 



2. The Lugo bolide 

On 1993 January 19 at 00:33:29 UT a large meteoroid 
entered the atmosphere approximately over the town of 
Lugo, in Emilia Romagna, Italy. The impact was recorded 
by the National Research Council (CNR) forward-scatter 
meteor radar and by six seismic stations, three belong- 
ing to the Microseismic Network of Ferrara (Pontisette, 
Ca Fornasina, Fiorile d'Albero) and the others to the Na- 
tional Institute of Geophysics (Barisano, Santa Sofia, Pog- 
gio Sodo). The event was also observed by several eyewit- 
nesses, as it lit an extremely large area (almost all of Italy), 
and they reported a visual magnitude in the range —22 to 
—25. Preliminary calculations were carried out based on 
the eyewitness reports, although they were fragmentary 



and sometimes contradictory (Cevolani et al. 1993, Kor- 
levic |1993 ). Only at a later time we found seismic data 
which enabled us to infer the location of the explosion 



(Cevolani et al. 1994). This analysis indicated that a me- 
teoroid of initial radius in the range 1.5 3 m impacted 
the Earth atmosphere at a velocity of about 26 km/s, with 
an inclination of the trajectory to the horizon of 8° 20°. 
By means of the seismic data, it was possible to calculate 
the height (30 ± 3 km), latitude (44°.48 ± 0°.01 N) and 
longitude (11°.91 ± 0°.01 E) of the explosion. 

3. The reanalysis: Aerodynamics 

Here, we will assume that the only reliable data are those 
recorded by the seismic stations, which in general are a 
very useful tool for understanding this kind of airburst 



(e.g. Bcn-Menahem 1975). Therefore, we assume as valid 
the height, latitude and longitude of the explosion only, 
i.e. those data calculated from seismic data (Cevolani et 



al. 1994) 
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The aerodynamics of large meteoroid/small asteroid 
impacts has been studied by several authors, sometimes 
with special reference to the 1908 Tunguska explosion (e.g. 
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Fig. 1. Seismic plot recorded at the Pontisette station. 
Time starts at 00:36:37.3 UT. Furthe r plots of this type 
can be found in Cevolani et al. (1994). 



In order to calculate the flight path angle, we have to 
solve two equations: 



dh 
'dt 

dO 
'dt 



y • sin I 

cos 9 



V 



R + h 



(2) 
(3) 



where g is the gravity acceleration [m/s^], R is the Earth's 
radius (we assume R — 6367 km, for about 45° latitude), 
and 9 is the flight path angle, measured from the horizon- 
tal. We assume that the meteoroid lift can be neglected. 
For the Tunguska cosmic body, Chyba et al. (1993) as- 
sumed a lift value of 10~^ and found that its influence on 
the results of these calculations is only about 1%. With 
all these assumptions, we obtain that the flight path an- 
gle during the final part of the atmospheric trajectory was 
9 — 5.0° ± 0.3°. Again, we have some disagreement with 
the previous results (8° 20°). This is probably due to 
the uncertainty of visual observations in these conditions; 
for such an event the surprise can reduce significantly the 
skills and reliability of eyewitnesses. 



Ceple cha an d McCrosky |197q, Cep lecha et al. |199S , Chyba 

|199^ ). 



et al. 1992, Hills and Goda 



1993, Lyne et al. 



Al- 
though the details may vary, there is a consensus that a 
30 km explosion height is typical for a carbonaceous chon- 
drite or a cometary body. In the theory of Hills and Goda 
( |1993| ) the height of first fragmentation is calculated com- 
paring the stagnation pressure in front of the meteoroid 
{Pmax = PoVe) the mechanical strength S of the cosmic 
body. We rearrange the formula to evaluate the meteoroid 
speed (Ve): 



Ve = W — exp 
V Po 



he 



(1) 



where po is the atmospheric density at the sea level 
[kg/m'^], ft-e is the height of first fragmentation [km] and 
H is the atmospheric scale height (about 8 km). For 
the strength, we assume S — 10^ Pa, that is an in- 
termediate value between those appropriate for carbona- 
ceous chondrites and for cometary bodies. We obtain 
Ve = 18 ± 3 km/s, that is a value much lower than that 
derived previously (about 26 km/s). 

Observing the seismic plots (e.g. Fig.|l|), we can con- 
clude that there was a single explosion (for a comparison 
with nuclear explosions, see Pierce et al. |1971| ). There is no 
evidence of multiple explosions, as it should occur during 
multiple fragmentation. Thus, for the Lugo bolide, Eq. (|l|) 
can be used by assuming that the first fragmentation cor- 
responded to the airburst. 



4. The reanalysis: Explosion energy 

To obtain an estimate of the explosion energy, we can 
use the relationship for the maximum velocity of displace- 
ment of the solid rocks, obtained from studies on un- 
derground nuclear explosions (Adushkin and Nemchinov 



19941) . We rearrange their equation in order to calculate 



the energy, when the distance and the displacement veloc- 
ity are known: 



E = k- D^- 



V \12/7 



V240y 



(4) 



where E is the explosion energy in kton of TNT; D is 
the distance of the sensor from explosion [km]; v is the 
displacement velocity [mm/s]. This formula is valid for 
D < 100 km: in our case, seismic stations were located at 
distances smaller than 70 km. The coupling coefficient k is 
introduced to take into account that, in order to produce 
rock displacements, an airburst is less effective than an 
underground nuclear explosion (at least by a factor 100). 
Moreover, there is a difference in the effective energy, be- 
cause the explosion of a meteoroid in the atmosphere does 
not involve nuclear fission, and this contributes about an- 
other factor 10. Finally, there is some increase of the wave 
amplitude with the height of burst up to 40 km (Pierce 
et al. 1971), which typically exceeds a factor 2; we as- 
sume a power increase by a factor 5. Overall, we estimate 
/fc = 100 • 10 ■ i = 200. 

5 

We have data from six seismic stations (for a com- 
ple te set of plots and other information, see Cevolani et 
al. 



, but transfer functions are available only for the 
three stations belonging to the Microseismic Network of 
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Fig. 2. Fourier analysis of the Pontisette seismic plot. 
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Fig. 3. Fourier analysis of the Ca Fornasina seismic plot. 

Ferrara. We have performed a Fourier analysis of the wave- 
form and found a peak at 1.4 Hz, for both Pontisette and 
Ca Fornasina, corresponding to the airburst (see Figs. || 
and We have not taken into account data from the 
Fiorile d'Albero station, because they show a strong back- 
ground noise overlapping the shock wave and preventing 
a reliable Fourier analysis. 

The transfer function has a nominal value of 175 
mV-s/mm for all stations and for frequencies greater than 
2 Hz. Below the cutoff frequency, the transfer function is 
drastically reduced, down to a value of 10 mV-s/mm for 
0.5 Hz. For a frequency of 1.4 Hz, we have a transduction 



factor of 52 mV-s/mm. The final results of our calcula- 
tions for the explosion energy from the seismic data with 
Eq. are shown in Table |l|. 

Table 1. Explosion energy calculated from the seismic 
data. 



Station 


D [km] 


V [fim/s] 


E [kton] 


Pontisette 


59 ±3 


41.040 ±0.002 


14 ±2 


Ca Fornasina 


63 ±3 


35.369 ± 0.002 


13 ±2 



We consider a mean value of 14 ± 2 kton, that is 
(5.9 ± 0.8) X 10^^ J. It is worth noting that we might have 
obtained more accurate values, but the saturation of the 
Barisano sensor introduced an error of 9% in the burst 
height calculations (Cevolani et al. 1994), which propa- 



gates to our results. On the other hand, had we not con- 
sidered the Barisano data, the available data would have 
been insufficient for a meaningful analysis. 

When a cometary body or a carbonaceous chondrite 
enters the atmosphere, almost all the kinetic energy is 
released in the explosion. Then we can calculate the me- 
teoroid mass, taking into account that during the path 
preceding the explosion the cosmic body undergoes a lim- 
ited mass loss: 
2E 

m= — = (4±l).10Mkg]. (5) 

In order to calculate the visual magnitude of the air- 
burst, we have to solve the equation: 

dm 

where r is the dimensionless coefficient for the meteor lu- 
minous efficiency. This coefficient mainly depends on the 
meteoroid speed and is quite uncertain (Ceplecha and Mc- 
Crosky 1976). Some authors think that for very bright 
bolides r ranges from 10 to 30% (Brown et al. 1996, Mc- 
Cord et al., 199?;). Others assume r va lues betwe en 1. 5 



(6) 



and 6.1% (Borovicka and Spurny 1996 , Ceplecha 1996). 
Here we assume r = 4.5%. 

Moreover, we assume that the meteoroid dissipated al- 
most all of its energy within a scale height. Then, solving 
Eq. (H) for the time during which the meteoroid exploded, 
we obtain t — 0.8 s. The corresponding value for the 
airburst luminosity is (5± 1) • 10"'^^ J/s. In order to express 
the luminosity in terms of absolute magnitude (i.e., the 
magnitude as observed at a 100 km distance), we can use 
the equation: 



M = -2.5 -(logioL- 2.63), 



(7) 



where we have rearranged the classical relationship in or- 
der to use the SI unit system. From Eq. (|^) we obtain 
M = — 22.7±0.5, a value consistent with visual observa- 
tions (—22 — 25). We stress the importance of the coeffi- 
cient t: ass uming a value of 10%, as suggested by McCord 
et al. (|l995|) , we would obtain M ~ -24. 
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5. Further results and discussion 



Table 2. Summary on the properties of the Lugo bolide. 



It is also interesting to check how the results are sensitive 
to the assumed value of the strength S. If we take S = 
10^ Pa, that is typical for cometary bodies, we end up 
with a cosmic body with a speed of about 6 km/s and an 
inclination of 2°. The mass would be about 3 x 10^ kg and 
the absolute visual magnitude —21. The airburst would 
have been 31 s long. These values appear unlikely. Note 
that a final velocity of 6 km/s is very close to 4 km/s, 
which Ceplecha (1994) indicated as necessary to have a 
meteorite fall. But for Lugo no meteorite was recovered. 

We can summarize some features of the Lugo bolide: it 
had a grazing trajectory in the atmosphere, it was proba- 
bly a carbonaceous chondrite, but it exploded at a height 
higher than usual and with a single airburst, without frag- 
mentations. The recent discovery by the NEAR probe of a 
carbonaceous asteroid (253 Mathilde) with a very low den- 
sity (about 1300 kg/m^) suggests the existence of porous 
bodies (i.e. bodies with internal cavities) among aster- 
oids (Yeomans et al. 1997 ). If we assume that the Lugo 
bolide was a porous carbonaceous chondrite, we have a 
body which was probably stronger than a cometary frag- 
ment, but which could explode at a higher altitude than 
those typical for stony objects, because of its porosity. It 
is very likely that porosity increases the burst efficiency: 
when ablation removes the surface of the body, cavities 
may appear which increase the aerobraking and generate 
a sudden deceleration. The kinetic energy then is rapidly 
transformed into heat, so that the body bursts within a 
scale height. This is consistent with a single explosion, 
without multiple fragmentation, as indicated by seismic 
plots (see Fig. p. 



6. Conclusions 

The Lugo bolide has been reanalysed by taking into ac- 
count only the data recorded by seismic stations. We sum- 
marize the main inferred properties of the bolide in the 
following Table |[ 

We are now carrying out calculations on the orbit and 
the dynamical evolution this bolide, whose results will be 
available soon. However, from the analysis described here 
it appears likely that the meteoroid was a porous carbona- 
ceous chondrite, somehow similar in constitution to the 
asteroid 253 Mathilde. The porosity would have increased 
the braking and as a consequence the airburst occurred at 
a height higher than for a compact carbonaceous chondrite 
object. 



Apparition time ( U 1 ) 


1993 01 19 00:33:29 ±1 s 


Latitude of airburst^ 


44.48° ±0.01° N 


Longitude of airburst'' 


11.91° ±0.01° E 


Airburst height" 


30 ± 3 km 


Explosion Energy 


14 ± 2 kton 


Mass 


(4± 1) • 10^ kg 


Abs. Visual Magnitude 


-22.7 ±0.5 


Velocity 


18 ± 3 km/s 


Inclination'' 


5.0° ±0.3° 


Path azimuth''''' 


146.5° ±0.5° 



Calculated in Cevolani et al. (1994) 
Over the horizon. 
Clockwise from North. 
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